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ABSTRACT 

A simple two-level numerical model using the quasi-geos trophic 

forecast equations is developed. Equations are linearized and 

friction is introduced in the surface layer. Solutions are obtained 

numerically by using the initial value approach. Two wind profiles, 

2 

U = -U tanh y/y and U = U sech y/y , are used and these are known 
to be unstable. For each wind profile the growth rate is determined 
as a function of the wave number. Some observed features of easterly 
waves are reproduced in the numerical solutions. 
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TABLE OF SYMBOLS AND ABBREVIATIONS 



A kinematic eddy viscosity 

e 

a inflow angle 

8 derivative of coriolis parameter at y = 0 

o 

C specific heat at constant pressure 

P 

f^ coriolis parameter at y = 0 

g gravity 

O) dp/dt 

Y 

R gas constant 

p density 

T average temperature from some standard atmosphere 

a (R^T/p^g) (5 t/Sz + g/C ) 

C 

z height 
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1. INTRODUCTION 



Synoptic scale disturbances in the tropical atmosphere have long 
been recognized (Palmer, 1951; Riehl, 1954). In general these dis- 
turbances possess a ”cold-core" structure with the colder air and most 
active weather located east of the wave (Yanai and Nitta, 1968b). 
Vertical motions associated with these waves are generally upward east 
of the wave axis and downward to the west (Yanai and Nitta, 1967). The 
source of energy for these waves is not fully known, however, Charney 
(1963) argues that when condensation is absent the large scale motions 
tend to be quasi-horizontal and quasi-nondivergent . If this argument 
is correct then the barotropic instability that arises as the result 
of horizontal shear in the Intertropical Convergence Zone could be the 
source of energy for the disturbances (Schminke, 1968). 

Kuo (1949) studied the stability characteristics of a barotropic 
zonal current, and showed that the barotropic zonal current is stable 
if the gradient of absolute vorticity has the same sign in the zonal 
current. Charney (1964) discussed the case where horizontal shear 
instability is produced in the equatorial convergence zone when the 
zone is located away from the equator and the converging air masses 
carry their angular momentum with them. 

The purpose of this paper is to devise and test a numerical model 
that is structured realistically, that is, the perturbation has a con- 
figuration and behavior similar to observed easterly waves. The main 
feature desired is to have the associated weather, and therefore the 
vertical motion, upwind of the trough axis. 
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2. BASIS OF THE MODEL 



The assumption that the source of energy for the wave mechanism is 
barotropic instability leads to the selection of a basic current that 
is independent of height. Bar oc linicity has been omitted but the 
atmosphere remains stratified. Charney (1963) used a scale analysis to 
show that the barotropic vorticity equation governs tropical motions in 
the absence of condensation. In another paper Charney (1969) discusses 
the role of the phase velocity in energy propagation in an easterly 
regime. Following Schminke (1968) the assumption is made that the 
easterly waves can be roughly described by a two-level quasi-geos trophic 
model. It is recognized that the quas i-geos trophic approximation is not 
very accurate close to the equator, however, the zone of dynamic 
instability is assumed to be associated with the Intertropica 1 Conver- 
gence Zone which is some distance from the equator so that the 
approximation is relatively accurate. 

The type of wind profile chosen will play a large role in determining 
the stability of the disturbance. Jacobs and Wiin-Nielsen (1966) inves- 
tigated barotropic instability and found that there are several unstable 
modes. Since a basic zonal current with no change in the vertical is 
used in this model, the wave with no vertical motion would predominate 
(Jacobs and Wiin-Nielsen, 1966). Easterly waves are observed to have 
vertical motion, so a method must be found to introduce it. This is 
accomplished by introducing friction in the surface boundary layer. 

Thus vertical motion, which depends upon the stratification of the 
atmosphere, is now forced into the model. 
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3. THE FORECAST EQUATIONS 



A simple two-level model is constructed by dividing the entire 
atmosphere into four layers of constant pressure differential, Ap/2 
(fig. 1), numbered 0 to 4 from top to bottom. Vertical motion is 
assumed to be zero at the 



top of the 
surface is 



0 



Ap 



2Ap 

Fig. 1. 

atmosphere , 



CD = 0 



0 






1 



CD 



2 






3 



CD = - pgW 



Two-level model used for prediction, 
while the vertical motion term at the earth 



CO = - Pgw^. (3.1) 

Charney and Eliassen (1949) used the Ekman theory to derive an 

expression for w^ which is ^^2 

= I ^ sin 2a (3.2) 

o 

where is the kinematic eddy viscosity and a, the inflow angle. The 
surface geostrophic vorticity can be approximated by 

C4=C3=V^Y3. (3.3) 

Following the development of Thompson (1961), but using different 
notation, begin with the quasi-geostrophic vorticity equation 

|^V^Y+ Ik xVf.V (V^Y) + 0. (3.4) 

Notice that when CD = 0, the above equation becomes the barotropic 
vorticity equation, which is sufficient to describe the main instability. 
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Apply this equation at levels 1 and 3 giving 

N 2 I 2 SY, 0^2 

^ V + Ik . Vfj-V(7 ^ ^ = 0, (3.5) 

5 2, 2 ^^3 K ■ ^2^ 

fj V »3 + |k k W 3 -V(V T 3 ) + e„ -5^ + £„— - 0. (3.6) 

Smaller layers could be used if it is desired to stop at some intermediate 
point, say the tropopause, rather than including the entire atmosphere. 

Next consider the quas i-geostrophic first law of thermodynamics in 
the form 



ll" • v(i) + f-6>= 0 , (3.7) 



where 



R^T par ^ g_1 

~ + c J' 

p g 



(3.7a) 



When this equation is applied at level 2 the result is 

(Tj - 73 ) + Ik X V — j— • VO'j-'fj) j— = 0 . (3.8) 



Define the following quantities 



Y = 
m 



Y = 
Tt 



fi + Y 3 



Yl - Y 3 



which implies that 



Y, = + Y^ 

i m i 



Yo = Y - Y^. 
3 m T 



(3.9) 

(3.10) 

(3.11) 

(3.12) 



Here Y,j, is proportional to the layer thickness, therefore is a nmeasiuire 
of its mean temperature. Using these definitions add (3.5) and (3.6) 
and divide the result by 2 , obtaining 



3 - + Ik X VY -V(V^Y ) + |k X VY-V(V^Y.„) + 

ot m ' m m' ’ x .rj*. 






(IB. 



0 , 



- £ 



o 2^p 



(3.13) 
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For the second forecast equation subtract (3.6) from (3.5) and eliminate 



CO 2 using (3.8) 




(3.14) 



where 




(3.14a) 



These are the prediction equations for the model. They can be 
linearized by separating the flow into an east-west current that varies 
only in y, and a small departure from this flow. In this event it is 
possible to treat various waves in x independently, hence the fields 
may be defined as follows 



where k is the x wave number. 

Substitute expressions for and into (3.13) and (3.14), 
separate the various sine and cosine terms, neglecting all products of 
the quantities A through D, and equating coefficients of the cosine kx 
terms gives 



= E(y) H- A(y,t) cos kx + B(y,t) sin kx, (3.15) 



“ C(y,t) cos kx + D(y,t) sin kx. 



(3.16) 




ay 



(3.17) 



For the sine terms the result is 







(3.18) 
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Repeating the procedure for (3.16) the final equations are 

S 2 ^ 2 \ , r^E SE ^ 2 . 2. a^E 1 

^ - ck - j = ^ ^ - 77 dJ 



by 



by" 



and 



- 3^ Dk + K (A - C), 

oy 



S rb^T) ^,2 ^ 2'\ , r^E _ 2 . 2. bE b^C ^ 

^ - Dk - D^ j = k C (k + ^ ) - ^ _ + 



(3.19) 



.3 



By 



c] + 3q Ck + K - k^^ (B - D), 



By 



By 



(3.20) 



where 



f g / A 



1/2 



(3.21) 



2RT o 

It is now necessary to derive the computational equation for vertical 
motion starting with 

2 r^'^T 

C0o=T^|-^+u ~^+v -^1, (3.22) 



where 



u 

m 



1 

f 



by 

m 

by ’ 



(3.23) 



V 

m 




f bx 



(3.24) 



o 

Substituting equations (3.23), (3.24), (3.15) and (3.16) into (3.22), 
collecting terms, and again neglecting all products of the quantities 
A through D results in 



“2 - ■'* Cl? - r If 0 

sl» kx + c)]. <3.25) 

o ^ 

This is the equation used in the model for the computation of vertical 
motion. The prediction equations and the vertical motion equation are 
now in the form to be used in the model. 
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4. BOUNDARY CONDITIONS AND FINITE DIFFERENCING SCHEME 



The finite differencing scheme used is illustrated below with a 
sample variable N; 



where i is the grid index and H is the distance between grid points. 

Centered time differences are used for all quantities except those 
involving friction. The frictional terms are computed at time (t - At). 
The first step in all cases is a forward time step. The second order 
equations for time tendencies are solved by the exact method of Richtmyer 
(1957, p. 101). 

Yanai and Nitta (1968a) studied the problem of finite difference 
approximations in solving dynamic instability problems of non-divergent 
barotropic currents. They showed that for both symmetric and anti- 
symmetric zonal currents the exact boundary conditions can be replaced 
by rigid boundary conditions at a distance equal to the half width of 
the shearing wind belt. At any boundary distance larger than this value 
the numerical and theoretical values tended to have a constant difference. 
Therefore rigid boundaries are placed at y = -w/2 and y = w/2 , where w is 
the total width over which computations are made. With the exception of 
E (the basic current) all variables have the boundary condition A = B = 
C=D=0aty= -w/2 and at y = w/2. When second derivatives are 
required on the boundary they are set equal to zero. 




(4.1) 






(4.3) 
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+ w/2 
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5. WIND PROFILES AND INITIAL CONDITIONS 



Schminke (1968) investigated and discussed the structure of easterly 

waves from a similar model using one wind profile (U = -UQtanh(y/y )) 

and one non-dimensional wave number (l/2y^). It is proposed to expand 

2 

this to include another profile (U = U^sech y/y^). Additionally the 
behavior of disturbances at various non-dimensional wave numbers under 
both profiles will be considered. 

Charney (1963) shows by scale analysis that in the absence of 
condensation, flow in the tropics is governed by the barotropic 
vorticity equation. The assumption is made that the source of energy 
for easterly waves must be barotropic instability. Since potential 
energy is assumed not to be an important energy source the basic 
temperature field may be set equal to a constant and is invariant. 
However, small temperature fluctuations are observed with easterly 
waves, so a two level model is needed (Schminke, 1968). 

If the waves are to grow, the zonal wind profile must be baro- 
tropically unstable. Following the suggestion of Charney (1963) that 
the shear layer near the Intertropica 1 Convergence Zone may be formed 
by transport of angular momentum, the unstable profile chosen to 
approximate this is 

U = -UQtanh y/y^, (5.1) 

illustrated in figure 2. This seems most appropriate during monsoonal 
conditions. The second profile chosen, 

U = U^sech^y/y^ , (5*2) 

is probably more applicable to mean flow conditions, and is illustrated 
in figure 3. 
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Garcia (1956) found that the profile given in (5.1) is barotropica lly 
unstable when 0 < k < Betchov and Criminale (1967) have analyti- 

cally, and by use of the eigenvalue approach, investigated this profile 
throughout the range mentioned above and found the maximum instability 
at a non-dimensional wave number of 0.45 (ky^ = 0.45). With these 
results in mind this paper is confined to the range 0.25 to 0.60 for 
the non-dimensional wave number when dealing with the hyperbolic 
tangent profile (5.1), Betchov and Criminale (1967) found a rather 
flat zone of maximum instability for the hyperbolic secant squared 
profile (5.2) in the range 0.80 to 1.05 for the non-dimensional wave 
number. Some computations were done in the lower range (0.25 to 0.60), 
but the majority of the experiments with this profile were confined to 
the range 0.70 to 1.15. Results of experiments with both profiles will 
be discussed in a later section. 

With the wind profile independent of height and with no friction Y 



m 



and should behave independently. Jacobs and Wiin-Nielsen (1966) 

showed that the growth rate for the Y field is greater than that for 

m 

the Y^ field. The thermal field Y^ is zero in the initial state but 

it changes as a result of the Ekman friction term. Y is introduced 

m 

into the model as a random disturbance in y. Introduction of surface 
friction leads to the coupling of the two fields. Thus at time t = 0 



C = 0 
D = 0 



or 



random disturbance 


(y), 


(5.3) 


random disturbance 


(y), 


(5.4) 


0 




(5.5) 


0 




(5.6) 


UoYolnCcoshCy/yo)) , 




(5.7) 


U y tanh y/y . 
0 0 *^ 0 




(5.8) 
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6. COMPUTATIONAL PROCEDURE 



A list of the constants used in the computations is given below; 

y = 200 kilometers, 
o 

w = 2000 or 4000 kilometers, 

U =10 meters per second, 

= 5 X 10 ^ per second, 

a =0.8 meters per second squared per centibar squared, 

13 = 0.0 or 2.29 x 10 per meter per second, 

Ap = 50 centibars, 

A^ = 10 meters per second squared, 

a =22,5 degrees , 

2 -12 

jj. = 2.49 X 10 per meter squared, 

H = 25 or 50 kilometers. 

At = 0.5 or 1.0 hours. 

When w = 2000 km and H = 50 km a forecast period of 14 days was used, 
in all other cases the forecast period was 28 days. A sample program 
as developed by Schminke (1968) and modified to fit the IBM 360 computer 
is included as Appendix 1 (note: this program is for w = 2000 km and 

H = 25 km). 

The random disturbance introduced into the ¥ field, (5.3) and 

m 

(5.4), was generated by a random number generation program, then read 
directly into the computer. The same disturbance field was used for 
all experiments with both variations of 3> H, At, and all variations 
of the non-dimensional wave number. After a few days the disturbance 
adjusted to an exponential growth. For convenience, growth rates were 
computed over the last three days of the forecast period. Once the 
field has adjusted it does not matter when computed, rates are always 
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the same. This was done by assuming the amplitudes could be written as 



= Sexp (U^n/y.^) 

X2 = Sexp (U^n/y^) t^. 



( 6 . 1 ) 

( 6 . 2 ) 



where n is the growth rate and S is a constant. By forming the ratio 



the growth rate can be shown to be 



n = 



y In (x /x ) 

O Z 1 

■ 



( 6 . 3 ) 
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7 . RESULTS 



A series of experiments was conducted under varied conditions'. In 
general the total distance over the grid was 2000 km, but one series was 
done over a total distance of 4000 km. Grid mesh was either 50 km or 25 
km with the smaller distance producing better results. When the grid 
mesh was 25 km a time step of 0.5 hours was used, however the majority 
of experiments were run at t = 1.0 hours. All cases were run with 
3=0, then with 3 ” 2.29 x 10 meters per second. Many combinations 
of friction, wave number and wind profile were run, and selected ones 
are illustrated in figures 4 through 19. 

Growth rates for the disturbances with both profiles and for various 
non-dimensional wave numbers are shown in figures 4 through 11. Betchov 
and Criminale (1967) found a maximum growth rate of approximately 0.19 
for a non-dimensional wave number of 0.45 using the hyperbolic tangent 
profile. With rigid boundaries, no friction and a fairly large grid 
size (H = 25 km) through the maximum shear zone (width approximately 
400 km), the result (fig. 7) was a growth rate of 0.165 at a wave 
number of 0.45. According to Yanai and Nitta (1968a) the best results 
should be obtained with a minimum of 20 grid steps across the zone of 
maximum shear. Therefore the results for the hyperbolic tangent pro- 
file are considered good, and the trend upward toward the Betchov and 
Criminale (1967) results is evident as the grid size is decreased 
(compare figures 5 and 7). The inclusion of the 3 parameter had no 
significant effect on the growth rates giving only a slight decrease 
of 0.01 or less. 
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The hyperbolic secant squared profile results were equally 
encouraging. Just as Betchov and Criminale (1967) lead one to expect, 
a relatively broad zone of maximum growth is found reaching from 0.85 
to 0.95 non-dimensional wave numbers. It is difficult to ascertain just 
what maximum rate Betchov and Criminale (1967) found, however, results 
of these experiments show rates that are significantly smaller than the 
rates for the hyperbolic tangent profile. The inevitable conclusion is 
that the hyperbolic secant squared profile is not quite as unstable as 
the hyperbolic tangent profile. As before, decreasing the grid size 
in the maximum shear area led to larger growth rates. Figures 9 and 10 
illustrate this very well. 

In figures 12 through 15 the amplitude and phase relationships of 

Y 5 and cn are depicted for the wave number of maximum growth, 
m T 

Amplitudes of Y , Y„ and CD show maxima in the region of greatest 
m T 

horizontal shear. The influence of the frictional term is clearly 
seen in the relative amplitudes of Y^ and Y^ • The three fields show 
close phase relationships and all phases change rapidly through the 
shear zone. It appears to be a safe assumption that the temperature 
and pressure fields are in close relation for both levels. The omega 
amplitude (fig. 14) shows a distinct minimum in the center of the 
field but sharp maximum peaks on each side in the maximum shear zone. 
Phase relations between Y^ and CD might almost be termed ideal with the 
sinking motion north of the shear zone about 70 degrees upwind of the 
amplitude ridge, placing the convergence in the area where it is 
normally observed in easterly waves. 

Addition of the 0 parameter made small insignificant changes in the 
amplitude; phase relationships maintained the same basic pattern as 
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without 3« Variation in the non-dimensional wave number over the range 
0.25 to 0.60 left the pattern basically unchanged while producing the 
growth rate changes discussed above. 

Figures 16 through 19 illustrate the behavior of the hyperbolic 
secant squared profile. Note that far exceeds the other components 
in amplitude (fig. 16) and that the maximum is symmetrical around the 
center of the maximum shear zone. fs very small indicating that 

friction is probably more effective with the hyperbolic secant squared 
profile. The lower field, Y^ > tends to lead the upper field in west- 
ward movement. This appears reasonable and is observed in tropical 
disturbances of this nature (Riehl, 1954, pp. 215). This gives the 
entire system a tilt toward the east with increasing height; this tilt 
is not noticeable with the hyperbolic tangent profile. 

The omega amplitude is very sharp peaked and symmetric about the 

center. Velocities, from about + y to the boundaries, are very small. 

■” o 

In this area the convergence appears to be about 100 degrees out of 
phase with the trough. Near y^ = 0, along the Intertropica 1 Conver- 
gence Zone, the convergence actually leads the trough. 

Use of a finite value for the 3 parameter made no significant 
difference in the results with this profile. 
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~l 1 1 1 1 J— 

0.30 0.40 O.iO 

Fig. 4. Growth rate for hyperbolic tangent profile, 
w = 2000 k"i, H = 50 km, with friction. 
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0.60 
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0.30 0.40 0.50 

Fig. 5. Growth rate for hyperbolic tangent profile, 
w = 2000 kni, H = 50 km, no friction. 



0.60 
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Fig. 6. Growth rate for hyperbolic tangent profile, 
w = 2000 km, H = 25 km, with friction 
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Fig. 7 
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0.40 



I 

0.50 



-r- 

0.60 



kv 
■ o 

Growth rate for hyperbolic tangent profile, 
w = 2000 km, H = 25 km, with no friction. 
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“ 1 — 

0.30 

Fig. 8 



T 



“T" 

0.40 



“I r 

0.50 



Growth rate for hyperbolic tangent profile, 
w = 4000 km, H = 50 km, with friction. 
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0.60 
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.110 



“I — 

0.80 

Fig. 9 



— I 1 1 1 1 1 r 

0.90 1.00 1.10 

Growth rate for hyperbolic secant squared profile, 
w =* 2000 km, H - 50 km, with friction. 
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Fig. 10. Growth rate for hyperbolic secant 

squared profile, w » 2000 km, H =■ 25 km, 
with friction. 




ky 

•^o 



Fig. 11. Growth rate for hyperbolic secant 
squared profile, w ■ 4000 km, 

H ■ 25 km, with friction. 
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Fig, 12, Amplitudes of Yp and Y^ for hyperbolic tangent 

profile plotted against arbitrary scale, w = 4000 km, 
H “ 50 km, and wave number * 0.45. 



3.0 




Fig. 13. Phase relationships corresponding to amplitudes of fig. 12. 
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Fig. 14. Amplitude of omega, arbitrary scale for same case as figs. 
12 and 13. 



32 



a> 




Fig. 15. Phase relationship between Y, and cu, same conditions 
as figs. 12, 13, and 14. 



33 




Fig. 16. Amplitude of , and Y^ with hyperbolic 

secant squared profile, w =• 2000 km, H = 25 km, 
ky^ ■ 0.95, and scale arbitrary. 
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Fig. 17. Phase relationship for , Y^ , Y_ for same conditions 
as fig. 16. 
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Fig. 19. Phase relationship between omega and for conditions 
in fig. 16. 
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8. CONCLUSIONS 



The easterly wave model developed appears to be reasonable if the 
assumption of barotropic instability as the energy source is accepted. 
Basic zonal flow has been introduced by two profiles, both independent 
of height. The temperature variations arose from the action of the 
Ekman friction terms. Waves developed by the model behaved in many 
respects as easterly waves have been observed to act. Results show 
that the growth rate of disturbances is a function of wave number, there- 
fore wave length, as well as other factors. Separation of the 
boundaries to greater distances made some small changes but the 
greater effect is obtained by decreasing the grid size, which also 
required changing the time step size, through the shear zone. 

Vertical motions produced by the hyperbolic tangent profile are 
very similar to those observed in easterly waves, in that the maximum 
sinking motion is upwind of the lower level ridge and the maximum 
rising motion is upwind of the lower level trough. The disturbance 
maintains its phase relationship in the vertical so that no tilt is 
observed. With the hyperbolic secant squared profile the amplitude 
of the vertical motion outside the maximum shear zone is small and 
again the rising motion is upwind of the trough and sinking upwind 
of the ridge. Within the area of large horizontal shear the con- 
vergent area moves to a position downwind of the low level trough. 

It is not known whether this is observa tiona lly verified along the 
Intertr opica 1 Convergence Zone. Using this profile, the hyperbolic 
secant squared, the model has developed a solution that indicates a 
definite eastward tilt with height. This tilt has been verified by 
observations. Perhaps some combination of the two profiles would be 

more realistic in depicting actual waves. 
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A more complete investigation of the effect of friction on the 
symmetric secant squared profile appears warranted. There is, also, 
a definite need to include the effects of condensation in the model. 

Non-linear effects, their role in changing the mean wind profile 
and the structure of the disturbance, are very important and must be 
considered. Use of the primitive equations or the balance equations 
would seem to offer the most effective means of incorporating these 
effects. They appear to be crucial if an easterly wave is to grow 
into a tropical storm. 

Always of utmost importance is the need for more detailed and 
extensive observations, but nowhere is this as true as of the tropics. 
No model can be developed and adequately tested without detailed 
observations . 
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THIS CUTS TIME STEP IN HALF FOR FIRST FORECAST TIME 
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